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SHOCKS AND CAVITIES FROM MULTIPLE OUTBURSTS IN THE GALAXY GROUP NGC 5813: A 

WINDOW TO AGN FEEDBACK 

S. W. RANDALl\ W. R. FORMAnS S. GIACINTUCCI^■^ P. E. J. NuLSEn\ M. SUN^ C. JONEsS E. CHURAZOV''•^ L. P. 

David\ R. Kraft^, M. Donahue^, E. L. Blanton'', A. Simionescu*, N. Werner* 

ABSTRACT 

We present results from new Chandra, GMRT, and SOAR observations of NGC 5813, the dominant 
central galaxy in a nearby galaxy group. The system shows three pairs of coUinear cavities at 1 kpc, 
8 kpc, and 20 kpc from the central source, from three distinct outbursts of the central AGN, which 
occurred 3 x 10^, 2 x 10'', and 9 x lO'' yr ago. The Ha and X-ray observations reveal filaments 
of cool gas that has been uplifted by the X-ray cavities. The inner two cavity pairs are filled with 
radio emitting plasma, and each pair is associated with an elliptical surface brightness edge, which we 
unambiguously identify as shocks (with measured temperature jumps) with Mach numbers of M 1.7 
and M « 1.5 for the inner and outer shocks, respectively. Such clear signatures from three distinct 
AGN outbursts in an otherwise dynamically relaxed system provide a unique opportunity to study 
AGN feedback and outburst history. The mean power of the two most recent outbursts differs by 
a factor of six, from 1.5-10x10'*^ erg s~^, indicating that the mean jet power changes significantly 
over long 10^ yr) timescales. The total energy output of the most recent outburst is also more 
than an order of magnitude less than the total energy of the previous outburst (1.5 x 10^^ erg versus 
4 x 10"''^ erg ), which may be a result of the lower mean power, or may indicate that the most recent 
outburst is ongoing. The outburst interval implied by both the shock and cavity ages (~ 10'' yr) 
indicates that, in this system, shock heating alone is sufficient to balance radiative cooling close to the 
central AGN, which is the relevant region for regulating feedback between the ICM and the central 
SMBH. 

Subject headings: galaxies: active — galaxies: clusters: general — galaxies: groups: individual 
(NGC5813) — galaxies: individual (NGC5813) — X-rays: galaxies 



1. INTRODUCTION 

A major result from early Chandra and XMM-Newton 
observations was that the amount of gas in cool core clus- 
ters that cools to low temperatures is less than predicted 
by classical radiative cooling models (David et al. 2001; 
Peterson et al. 2001; Peterson & Fabian 2006). The im- 
plication is that the central gas must be re-hcated. The 
source of this heating, and understanding when and how 
it takes place, has recently been a major topic of study. A 
promising candidate is feedback from energy injection by 
the central AGN of the cD galaxy (e.g., Churazov et al. 
2001; Churazov et al. 2002; for a review see McNamara 
& Nulsen 2007). However, the details of this interaction, 
and how the energy is transferred from the jcrts to the 
ambient ICM, are poorly understood. Galaxy groups 
provide an excellent opportunity to study heating and 
other non-gravitational processes in the ICM. Although 
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not as X-ray luminous as clusters, the effects of heating 
are more readily apparent in groups, due to their lower 
gas temperatures, masses, and central densities. For ex- 
ample, the gas to total mass fraction in groups ranges 
from 0.02-0.07, with a scatter of ~ 2 at a fixed tem- 
perature within r25oo- The scatter is tightly correlated 
with the central entropy (Gastaldello et al. 2007; Sun et 
al. 2009), reflecting the greater role of non-gravitational 
processes in the centers of groups as compared to clus- 
ters. 

In this paper, we report on Chandra observations of 
NGC 5813 (UGC 09655), a bright (My = -22.01, Lauer 
et al. 2007) El galaxy. It is the central dominant mem- 
ber of a subgroup (which we shall call the NGC 5813 
galaxy group) in the NGC 5846 galaxy group (Mahdavi 
et al. 2005). NGC 5846 and NGC 5813 have a projected 
separation of 79.7' (740 kpc). This group is relatively iso- 
lated, lying well off the plane of the Local Supercluster. 
Both the NGC 5813 and NGC 5846 galaxy groups are 
members of the ROSAT-ESO Flux-Limited X-ray (RE- 
FLEX) galaxy cluster catalog (Bohringer et al. 2004). 
Detailed Hubble Space Telescope (HST) observations re- 
veal that NGC 5813 contains a relatively undisturbed 
dusty circumnuclear disk (Tran et al. 2001), suggesting 
that this galaxy has not recently experienced a major 
merger. Emsellem et al. (2007) classify NGC 5813 as a 
"slow rotator" galaxy which, they argue, represents the 
extreme evolutionary end point reached in deep poten- 
tial wells, consistent with this object being a dynam- 
ically old galaxy at the center of a galaxy subgroup. 
NGC 5813 contains a 2.8 x IO^Mq supermassive black 
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hole (SMBH), and an associated AGN that is a known 
radio source (e.g., Balmaverde & Capetti 2006; Maghoc- 
chetti & Briiggen 2007). 

We report here on a 150 ksec combined observation of 
NGC 5813 with the Chandra X-ray Observatory, an anal- 
ysis of archival multi-frequency Very Large Array ( VLA) 
observations, Ha observations with the Southern Astro- 
physics Research Telescope [SOAR), and on some initial 
results from low frequency Giant Metrewave Radio Tele- 
scope (GMRT) radio observations. We focus on three 
main results: 

1. The ICM in NGC 5813 shows clear signatures from 
three distinct AGN outbursts, with three pairs of 
roughly coUinear cavities and unambiguous shocks 
with measured temperature jumps associated with 
the inner and intermediate cavities. 

2. The mean power of the two most recent outbursts 
differs by a factor of six, from 1.5-10x10*^ erg s~^, 
even in this otherwise dynamically relaxed system, 
indicating that the mean jet power varies over long 

10^ yr) timescales. 

3. The heating from shocks alone is sufficient to offset 
radiative cooling locally (within at least 10 kpc), 
without requiring the internal energy of the X- 
ray cavities. The heating is roughly isotropic, and 
strongest near the AGN where the shock Mach 
numbers are larger, which is the region of inter- 
est for regulating feedback between the ICM and 
the central SMBH. 

We assume an angular diameter distance to NGC 5813 
of 32.2 Mpc (Tonry et al. 2001), which gives a scale 
of 0.15 kpc/" for flo = 0.3, flf, = 0.7, and Hq = 
70 km Mpc~^. All uncertainty ranges are 68% con- 
fidence intervals (i.e., Icr), unless otherwise stated. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Chandra Observations 

NGC 5813 was originally observed with Chandra on 
April 2, 2005, for 49 ksec (ObsID 5907) with the Chan- 
dra CCD Imaging Spectrometer (ACIS), pointed such 
that the galaxy was centered on the back-side illuminated 
ACIS-S3 CCD (the ACIS-Sl as weU as the front-side 
illuminated ACIS-I3 and ACIS-S2 CCDs were also ac- 
tive). It was subsequently observed for 100 ksec on June 
5, 2008 (ObsID 9517) with the same chip configuration. 
These data were reduced using the method described in 
Randall et al. (2008). All data were reprocessed from 
the level 1 event files using the latest calibration files 
(as of CIA04.2). CTI and time-dependent gain cor- 
rections were applied. lc_CLEAN was used to remove 
background flare^ The mean event rate was calculated 
from a source free region using time bins within 3a of 
the overall mean, and bins outside a factor of 1.2 of this 
mean were discarded. There were no periods of strong 
background flares. The resulting cleaned exposure times 
were 48.7 and 99.6 ksec, respectively. 

Diffuse emission from NGC 5813 fills the image FOV 
for each observation. We therefore used the CALDEF°] 



blank sky background files appropriate for each obser- 
vation (including the new "period E" files for the more 
recent observation), normalized to our observations in 
the 10-12 keV energy band. To generate exposure maps, 
we used a MEKAL model with kT = 0.7 keV, Galac- 
tic absorption, and abundance of 30% solar at a redshift 
z = 0.006578, which is consistent with typical results 
for the extended emission from detailed spectral fits (see 
§11. 

The exposure corrected, background subtracted, 0.3- 
2 keV Chandra image is shown in Figure [T] To enhance 
the visibility of the diffuse emission, bright point sources 
were removed, and the regions containing point sources 
were "filled in" using a Poisson distribution whose mean 
was equal to that of a local annular background region. 
A close-up of the core in Figure [2] shows small-scale struc- 
ture in the center, while the more heavily smoothed im- 
age in Figure [3] shows structure in the fainter outer re- 
gions. We discuss the main features in these images in 

^11 X-ray spectra were fitted in the 0.6-3.0 keV band 
and grouped to a minimum of 40 counts per spectral 
bin. The absorption was fixed to the Galactic value of 
Nh = 4.37 X 10^° cm-2 (Kalberla et al. 2005). Vary- 
ing the absorption from the Galactic value did not sig- 
nificantly improve the fit. Anders & Grevesse (1989) 
abundance ratios were assumed throughout, unless oth- 
erwise stated. Temperature maps were derived using the 
method of Randall et al. (2008). For each temperature 
map pixel, we extracted a spectrum from a circular re- 
gion containing a minimum number of net counts (af- 
ter subtracting the blank sky background) in the 0.6 - 
3.0 keV band. The resulting spectrum was fitted with an 
absorbed APEC model using xspec, with the abundance 
allowed to vary. 

2.2. Radio Observations 

NGC 5813 was observed with the CMRT aX 235 MHz, 
as part of a larger, ongoing project (Giacintucci et al. 
2009; Giacintucci et al. 2010). The observations were 
carried out in August 2008 for a total of 100 minutes on 
source. The NRAO Astronomical Image Processing Sys- 
tem (AIPS) package was used for the data reduction and 
analysis. The visibilities were inspected and edited to 
identify and remove bad data. After the initial calibra- 
tion, phase-only self-calibration was applied to remove 
residual phase variations. Multi-field imaging was im- 
plemented in each step of the data reduction. The rms 
noise level (Icr) achieved in the image is 0.3 mJy/beam. 
We refer to Giacintucci et al. (2010) for a detailed de- 
scription of the data reduction and analysis. 

We also include results from an analysis of archival 
VLA observations at multiple frequencies. Data calibra- 
tion and imaging were carried out in AIPS following the 
standard procedure (Fourier-transform, Clean and Re- 
store) . Phase-only self-calibration was applied to remove 
residual phase variations and improve the quality of the 
image. The properties of the VLA and GMRT radio ob- 
servations are summarized in Table [TJ 

2.3. Ha Observations 
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Narrow-band Ha imaging observations were made with 
the SOAR Optical Imager (SOI) on July 6, 2008 (UT). 
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The night was clear and photometric. Two CTIO 
narrow-band filters were used, 660075-4 for the Ha-|-[NII] 
lines and 6129/140 for the continuum. Three 15-minute 
exposures were taken with the 660075-4 filter and three 
12-minute exposures were taken with the 6129/140 filter. 
Spectroscopic standard stars were EG274 and LTT7987. 
More detail on the SOI data reduction can be found in 
Sun et al. (2007). 

3. X-RAY CAVITIES AND SHOCKS 

As we discuss in detail in this paper, the X-ray im- 
ages show signatures from three episodes of AGN ac- 
tivity. There are two clear pairs of cavities distributed 
coUinearly (in a direction roughly parallel to the minor 
axis of the elliptical optical isophotes) and symmetri- 
cally about the galaxy center, each pair associated with 
a shock. In addition, we argue for a third pair of cavi- 
ties at larger radii (along roughly the same line) and an 
associated third surface brightness edge, which may be a 
weak shock. In detail, we see the following: 



1. A pair of inner cavities at 



the most recent AGN outburst about 3 x 10^ 



1.4 kpc, inflated by 
yr 

ago (Figure l2|. The cavities are surrounded by 
bright rims of emission, similar to what is seen in 
other systems (e.g., the Perseus Cluster Fabian et 
al. 2003, M87 Forman et al. 2007; M84 Finoguenov 
et al. 2008, Abell 2052 Blanton et al. 2009; see 
Figure^, though here we find the rims to be hotter 
than the ambient gas, whereas in other systems 
they are typically cooler. We identify the sharp 
edge in the rims, 1.4 kpc southeast of the AGN, 
as a shock, which we refer to as the 1.5 kpc shock. 



with a Mach number of M = 1.7 (see § |5.1[ ). The 
rims overlap to form an indented structure to the 
northwest. The northeastern inner cavity has an 
irregular morphology, possibly due to the wall of 
the cavity being "punched through" by the AGN 
jet. 

A pair of intermediate cavities at ^ 8 kpc, inflated 
by a previous AGN outburst about 2 x 10^ yr ago 
(Figure[l] Fi gurepl). The southwestern cavity has a 
regular morphology, while the northeastern cavity 
is more extended in the radial direction and may 
be a double "Russian doll" cavity (e.g., as in M84, 
Finoguenov et al. 2008). The intermediate cavi- 
ties lie just inside a sharp, elliptical edge, which 
we identify as a shock (and refer to as the 10 kpc 



shock) with M = 1.5 (sec § 5.1) 



A faint pair of outer cavities, at ^ 20 kpc, inflated 
by aprevious outburst about 9 x 10^ yr ago (Fig- 
ure ^ Figure E| . The northeastern outer cavity 
is surrounded by a rim of brighter emission, while 
the so uthw estern outer cavity is a weaker feature 
(see § 5.3). The outer cavities are associated with 
a weak outer edge-like feature at ^ 25 kpc (Fig- 
ure [sj, just past the outer edges of the cavities. 
Thisieature may represent an old shock, or a tran- 
sition region from the galaxy at mosp here to the 



extended group atmosphere (see § 5.2) 



X-ray luminosity of Lx = 1.6 x 10^^ erg s~^ in the 
0.3 ~ 12.0 keV rest frame energy band. Although 
this source is quite faint, we identify it as the AGN 
that has inflated the X-ray cavities and driven the 
shocks in the ICM due to its central location and 
detection in the radio (see 



4.2.5) 



4. A central point source, offset ~ 0.5 kpc southeast of 
the axis defined by the cavities (Figure [2]), with an 



Since it is difficult to show all the cavities simultaneously 
in the same image, we divided the X-ray image by a /?- 
model to better show surface brightness fluctuations over 
a wider dynamic range. The resulting image is shown in 
Figure [4j with the cavities indicated by overlaid regions 
for clarity. The overall morphology of NGC 5813 is re- 
markably symmetric and regular, suggesting that AGN 
feedback maintains a near "steady state" through regular 
outbursts in an otherwise undisturbed system (consistent 
with optical studies that also conclude NGC 5813 is dy- 
namically old, Tran et al. 2001; Emsellem et al. 2007). 

3.1. Radio Emission from the X-ray Cavities 

Radio contours from our 235 MHz GMiJT observations 
(green) and archival 1.36 GHz VLA observations (blue) 
are shown overlaid on the Chandra image of the core 
in Figure [5] Extended radio emission at 1.36 GHz fills 
the inner cavities, while the 235 MHz emission overfills 
the cavities, extending along the axis of symmetry of 
the cavities (this extension is real, and not due to the 
larger beam size at 235 MHz). The intermediate cavi- 
ties are filled with 235 MHz radio emission, but show no 
emission at 1.36 GHz. The intermediate cavities are also 
not detected at 1.4 GHz in the NRAO VLA Sky Survey 
{NVSS, which has a beam size of 45") or at 1.36 GHz in 
archival VLA C array configuration observations (see Ta- 
ble [T]). The outer cavities, which are outside the FOV of 
Figure [5j do not show any detected radio emission. The 
radio images are therefore qualitatively consistent with 
what is expected for intermittent AGN outbursts, where 
the electrons contained in the cavities age due to syn- 
chrotron, inverse Compton, and adiabatic losses as the 
cavities rise buoyantly after the outburst phase. The cen- 
tral cavities contain recently accelerated electrons, which 
emit at high and low frequencies, while older cavities 
contain older electron populations with fewer energetic 
particles and weaker high frequency emission. 

4. THE THERMAL STRUCTURE OF THE GAS 
4.1. Temperature Map 

The X-ray image (Figure Il| shows complicated struc- 
ture in the ICM, which fiUs the FOV. To study the ther- 
mal structure of the ICM, we generated a temperature 
map, requiring 1500 net counts per extraction region. 
The resulting temperature map, with the X-ray cavity 
regions overlaid, is shown in Figure[6j The corresponding 
pseudo-pressure and pseudo-entropy maps are shown in 
Figure [tI The extraction radii range from 2.8" (0.4 kpc) 
in bright regions near the core to 59" (9 kpc) in faint 
outer regions. The temperature uncertainties are be- 
tween 2% - 3% across the map. 

The temperature map shows that even in the pro- 
jected, effectively smoothed map, the hot (0.7-0.75 keV) 
10 kpc shocks are visible at the location of the promi- 
nent surface brightness edges. The pseudo-pressure map 
also shows large jumps across the edges, consistent with 
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these features being shock fronts. There is a trail of cool 
0.55 keV gas though the galaxy center, along the line 
defined by the X-ray cavities indicated in Figure [l] ter- 
minating at the edges of the intermediate cavities (we 
discuss this feature further in § 5.4 1. The kT ^ 0.65 keV 
gas extends to larger radii (out to ^ 27 kpc) in the east- 
northeast, coincident with the extension of diffuse emis- 
sion across the outer edge in Figure [3j East of this exten- 
sion, the temperature rises rapidly from about 0.65 keV 
to 0.75 keV over ^ 7 kpc. 

To study the detailed structure in the core, we made a 
higher angular resolution temperature map of the central 
region of NGC 5813. In addition to the finer spatial bin- 
ning, each extracted spectrum had only 1000 net counts, 
giving smaller extraction radii (and thus less smoothing 
in the map), ranging from 0.4 kpc to 1.8 kpc. The result- 
ing temperature map is shown in Figure |8j and the corre- 
sponding pressure and entropy maps in Figure[9j with the 
0.3-2.0 keV X-ray surface brightness contours overlaid. 
The bright rims surrounding the innermost bubbles are 
revealed to contain relatively hot {kT « 0.7 keV), high 
pressure gas, in contrast to the cool bubble rims seen in 
some other systems (e.g., the Perseus Cluster, Fabian et 
al. 2003; M87, Forman et al. 2007; Abell 2052, Blanton 
et al. 2009). The gas in the rims has likely been shock 
heated by a recent AGN outburst, which has rapidly in- 
flated the innermost cavities. 

4.2. Detailed Spectral Analysis 
4.2.1. Aztmuthally Averaged Profiles 

We produced projected radial profiles by fitting spec- 
tra extracted from concentric annuli, centered on the cen- 
troid of the diffuse emission at larger radii. Each annular 
bin was fitted with an absorbed APEC model, with the 
abundance allowed to vary. The resulting temperature 
profile is shown in the top panel of Figure |10| There 
is a temperature increase of ^ 0.05 keV at ~10 kpc, 
at the location of the surface brightness edges indicated 
in Figure [T] unambiguously identifying these features as 
shocks. Additionally, there is an inner temperature jump 
of ~ 0.06 keV at ^2 kpc, at the location of the bright 
rims around the inner bubbles, consistent with the core 
temperature map shown in Figure IS] Although the point 
source presumed to be the central AGN has been ex- 
cluded, the temperature profile shows a strong central 
spike. This is because the profile center is in the region 
of the hot overlapping rims of the central cavities, shown 
in Figure [2] 

4.2.2. Deprojection Analysis 

To determine the 3D structure of the ICM we per- 
for med a d eprojection analysis using concentric annuli as 
in § |4.2.1[ but with bins 2-3 times larger to provide ad- 
equate statistics for the deprojection analysis. We used 
the "onion peeling" method (employed, e.g., by Fabian 
et al. 1981, Blanton et al. 2003) to derive deprojected 
profiles. First, the projected temperature, abundance, 
and XSPEC normalization are determined by fitting an 
absorbed APEC model to the outermost annulus. Fits to 
spectra from annuli at smaller radii are then determined 
by adding an additional component for each outer annu- 
lus, with fixed temperature and abundance, and a nor- 
malization scaled to project from the outer to the inner 



annulus, assuming spherical symmetry. We note that this 
procedure does not correctly account for the uncertain- 
ties, since the contributions from outer shells are fixed. 
This method was adopted instead of simultaneously fit- 
ting spectra in all of the bins (e.g., with the projct 
model in xspec) to prevent spectra at small radii, where 
the spherical symmetry approximation is less accurate, 
influencing the fitting results at large radii. A compari- 
son of the two methods shows that while the best-fitting 
temperature profile is not significantly affected, the den- 
sity profile determined with the projct model differs 
somewhat from the profile we present here and shows a 
^50% increase in density with radius between 6-9 kpc 
(in the region of the intermediate bubbles and 10 kpc 
shock front). Thus, while the azimuthally averaged de- 
projected temperature profile is not significantly affected 
by the assumption of spherical symmetry, the density 
profile is sensitive to this assumption and is therefore 
only approximate. 

As a further check on the deprojection, we performed 
an independent analysis, with the data analyzed as de- 
scribed in Vikhlinin et al. (2005) and the deprojection 
method described in Churazov et al. (2008, 2010). This 
method accounts for emission at large radii by assum- 
ing a power law density profile for emission outside of 
the FOV and fitting the normalization of this component 
along with the fluxes from spherical shells. The derived 
temperatures agreed within the Icr confldence range, and 
the temperature jumps were fully consistent within the 
uncertainties. 

In the case of NGC 5813, measuring the deprojected 
abundance proflle is extremely difflcult. This is because 
at these relatively low gas temperatures (0.6-0.7 keV) 
the spectral resolution of the ACIS instrument is such 
that there is a strong degeneracy between line emission 
(which determines the abundance) and continuum emis- 
sion (which determines the electron density). When this 
effect is accumulated across several shells and projected 
onto the inner shells, the inner abundances are essentially 
indeterminate, leading to a highly uncertain density pro- 
file. Furthermore, determining abundances in multiphase 
gas is a known problem (e.g., Buote 1999; Rasia et al. 
2008). We therefore fixed the abundance at 50% solar, 
which is an average value from the projected proflle fits. 

The resulting deprojected (red triangles) and projected 
(black circles) temperature profiles are shown in the top 
panel of Figure |10[ alongside the corresponding density, 
pressure, and entropy profiles. The entropy is taken to be 

2 /3 

S = kTrio and the pressure P = nkT, where rie is the 
electron density and n = l.Sric- The most significant ef- 
fect of the deprojection on the temperature profile is the 
lower temperature of the gas at ~5 kpc (0.55 keV ver- 
sus 0.61 keV) once the projection effects from the outer 
shock-heated gas have been removed, as well as a signif- 
icant increase in the central temperature (from 0.65 keV 
to 0.7 keV). Note that in these plots the shock front 
edges are smeared over multiple bins since they are not 
spherically symmetric with respect to the galaxy center. 

4.2.3. Temperature Profile Across the Shocks 



rofiles (Fig- 

10) and projected temperature (Figure k3 I map show 
temperature rises at the location of the surface bright- 



The azimuthally averaged temperature 
ure 
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ness edges, characteristic of shock fronts. We therefore 
extracted spectra in sectors across these edges to bet- 
ter characterize the temperature jumps. Each sector was 
centered on the center of curvature defined by the cor- 
responding edge, which is not coincident with the cen- 
tral position used to extract the azimuthally averaged 
profiles. The extraction regions were truncated at small 
radii to avoid the complex structure in the central re- 
gion (additionally, the assumption of spherical symmetry 
breaks down at small radii since the centers of curvature 
are not coincident with the overall centroid of the diffuse 
emission). The projected and deprojected temperature 
profiles are shown in Figure [TT] The northwestern region 
shows higher overall temperatures and a larger temper- 
ature jump of ~ 0.15 keV as compared to ~ 0.1 keV in 
the southeast (where we estimate the size of the jumps 
by extrapolating the base-line deprojected temperature 
increase on either side of the temperature peak, roughly 
from r < 10 kpc and r > 18 kpc). The lower tempera- 
tures in the southeast are likely due in part to the east- 
northe ast e xtension of cool gas seen in the temperature 
map (§ 4.1), which overlaps with the extraction annuli. 



4.2.4. Total Diffuse Emission 

To accurately determine the weighted average prop- 
erties of the diffuse gas, we extracted and fitted spec- 
tra for the total diffuse emission within 2.9' (27 kpc). 
We initially fitted the spectra with a single APEC model. 
The resulting model showed residuals near the 1.8 keV 
Si and 2.46 keV S lines, possibly indicating non-solar 
abundance ratios. We therefore fit the spectra with a 
VAPEC model with the abundances of O, Ne, Mg, Si, S, 
and Fe allowed to vary independently (other elements 
were not well constrained and were fixed at 1/2 the 
solar abundance, except He which was fixed at solar). 
The resulting fitted abundance values were 0.8 - 1.0 so- 
lar, with Si and S having somewhat larger values than 
Ne and Mg, except for O (Zq = 0.35 ± 0.03) and Fe 
{Zpe — 0.65±0.03). Since the total diffuse emission is ex- 
pected to include emission from gas at multiple tempera- 
tures, we added a second vapec component to the model 
and tied individual abundances for the two components 
together. This model provided a much improved fit, with 
abundance values of Zq = 0.13lE] J]^, Znc = 0.59t;] J^^, 
ZMg = 0.65 ± 0.04, Zs; = 0.78 ± 0.04, Zs = 0.94 ± 0.08, 
and Zpe = 0.53 ± 0.03. The fitted temperatures were 
kTi = 0.351:^:;^? and fcTa = 0.67ll:^:;^j7. We tried includ- 
ing a power-law component in each of the above models, 
but in no case did this addition improve the fit (even if 
the photon index was fixed at a typical value of 1.5 for 
an unresolved LMXBs population) or tightly constrain 
the photon index. This suggests that emission from un- 
resolved LMXBs is not significant in this energy band. 
We note that the above results depend somewhat on the 
adopted abundance table. If we adopt the updated abun- 
dance table of Grevesse & Sauval (1998), we find that Fe 
is no longer under-abundant compared to the other ele- 
ments, with Zpe, Zs, and Zsi all roughly 0.7 solar, and 
Z]\fe and Zmq about 0.6 solar. O is still found to be 
under-abundant, with Zq = 0.16lo;o3- This difference 
in the fitted value for Zpe was noted previously, e.g., by 
Humphrey et al. (2004). 

As a further check of this result, we examined data 



from XMM-Newton RGS observations, which has bet- 
ter spectral resolution in the region of the O lines. The 
processing of the RGS data is described in Werner et 
al. (2009). We fitted the spectrum from a 1' wide re- 
gion centered on the core of NGC 5813 between 10-25 A 
with a collisionally ionized plasma (cie) model with two 
cooling flow components, one to account for gas cooling 
down to 0.4 keV and a second to account for gas cool- 
ing to lower temperatures. The abundances were con- 
strained to be equal across each model component. We 
found best-fitting abundance values of Zp^ = 0.60±0.06, 
Zo = 0.44 ± 0.06, and ZMe = 0.42 ± 0.06. Although this 
gives a higher Zo/Zp^ ratio (0.7 versus 0.2 from Chan- 
dra), O is still found to be under-abundant as compared 
to Fe. We conclude that there is evidence for non-solar 
abundance ratios in the diffuse ICM, with the most ro- 
bust result being a decreased O abundance relative to 
solar. Similar sub-solar values for Zq /Zp^ have been re- 
ported for other galaxy groups, massive elliptical galax- 
ies, and clusters of galaxies, and imply a greater rela- 
tive enrichment from Type la supernovae as compared 
to Type II (e.g., Finoguenov et al. 2000; Finoguenov et 
al. 2001). However, sub-solar ZojZp^ values are diffi- 
cult to reconcile with the larger Zsi/Zpf. and ZgjZp^ 
abundance ratios when one tries to apply SN enrich- 
ment models to explain the observed ICM abundances 
(see Humphrey et al. 2004 and references therein for a 
discussion). 

4.2.5. The Central Source 

NGC 5813 contains a central X-ray point source, vis- 
ible in Figure [2] We extracted a spectrum for the cen- 
tral source using an aperture with a radius of 1.5"and 
a background determined from a local annular region. 
The spectrum, which had 550 counts in the 0.6-5.0 keV 
band, was fitted with an absorbed power-law, giving a 
best-fitting photon index of F = 1.5^Q g. The nuclear 
X-ray luminosity is Lx = 1.6^05 x 10^^ erg in 
the 0.3 - 12.0 keV rest frame energy band (90% con- 
fidence ranges). We calculated the radio luminosity be- 
tween 10 MHz and 100 GHz using archival VLA A array 
configuration observations of the core at 1.49 GHz and 
4.86 GHz. This gave a spectral index of — 0.35 and 
a luminosity of Lfi = 1 x 10^^ erg s~^. Even though the 
spectral index indicates a relatively fiat spectrum for a 
core AGN, giving a larger estimate of the broad band lu- 
minosity than would a steeper spectrum, the calculated 
luminosity is still two orders of magnitude fainter than 
the faintest source in the sample of systems containing 
X-ray cavities given in Birzan et al. (2008). However, 
they mainly consider sources at the cores of galaxy clus- 
ters, which are higher mass systems than the subgroup 
we consider here. 

Although the central source is X-ray faint enough to 
be classified as an ultraluminous X-ray source (ULX) , we 
identify it as the AGN that has inflated the X-ray cavi- 
ties and driven the shocks in the ICM due to its central 
location and detection in the radio (see Figure [s]). Its 
low luminosity identifies this source as a low lummosity 
AGN, which typically have Lx < 1 x 10''^ erg s~^ (Ptak 
2001). This suggests that the source is either faint yet 
mechanically powerful, in a quiescent state after having 
recently had an outburst that inflated the inner X-ray 
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cavities, or heavily obscured (allowing the absorption to 
vary gave a value that was 4 times Galactic, although the 
fit was not improved and the absorption was consistent 
with the Galactic value to within la). 

5. DISCUSSION 

5.1. Structure of the Shock Fronts 

The hard band image, temperature and pressure maps, 
and deprojected temperature and pressure profiles iden- 
tify the prominent surface brightness edges around the 
intermediate cavities as shock fronts. To quantitatively 
study the structure of these shocks, we extracted the 0.3- 
2.0 keV surface brightness profiles from the sectors used 
to derive the temperature profiles shown in Figure 11 
Note that the profiles are defined by the centers of cur 
vature of the 10 kpc shocks, which do not coincide with 
the location of the central AGN. The profiles were then 
converted to integrated emission measure (lEM) profiles, 
using the temperatures from the projected temperature 
profile shown in Figure 11 with the abundance fixed at 
50% solar (roughly the av erage from projected profile fits 



in this region, se e §[4.2.2 1. The resulting lEM profiles are 
shown in Figure |l2| Each profile shows a sharp edge at 
~13 kpc. Following our previous work (Vikhlinin et al. 
2001; Randall et al. 2009a, 2009b), we fit the profiles 
by projecting a 3 dimensional density profile consisting 
of two power laws, connected by a discontinuous break, 
or "jump" . The free parameters were the normalization, 
the inner (a) and outer {/3) power law slopes, the position 
of the density discontinuity (fobroak); and the amplitude 
of the jump (A). The best- fitting model is shown as 
the solid lines in Figure [T2] with the fitted break radii 
indicated by the vertical oashed lines. The best fitting 
inner density jumps for the northwest and southeast sec- 
tors are A^^ = l.Tbtom and A^o = 1.69 ± 0.07. Us- 
ing the Rankine-Hugoniot shock jump conditions for a 
7 = 5/3 gas implies Mach numbers of M^w = 1.53 and 
Msc = 1.48, and temperature jumps by the same factor 
as the Mach numbers. The predicted temperature jumps 
are greater than the temperature jumps of ~ 1.2 detected 



in the deprojected temperature profiles (Figure 111. We 
discuss this discrepancy below. 

The core temperature map (Figure [8]) an d th e az- 
imuthally averaged temperature profile (r .gure 10 1 show 



evidence for shock heated gas in the bright rims sur- 
rounding the innermost bubbles, about 1 — 2 kpc from 
the central AGN. The X-ray image of the core (Figure [2| 
shows a sharp contrast between the bright bubble rims 
and the surrounding gas, with the sharpest edge to the 
southeast. We extracted the surface brightness profile 
in a 52° wide sector across the southeastern edge, cen- 
tered on the AGN, and fit the JEM profile with a broken 
power law density model as above. To the northwest, 
the morphology is more complicated, with the two rims 
meeting to form an indented structure. Since this struc- 
ture is not well-modeled by our assumption of spherical 
symmetry we did not fit the lEM profile to the north- 
west. There were too few counts to accurately measure 
the temperature profile across the southeast, so we as- 
sumed an isothermal gas with kT = 0.6 keV and 50% so- 
lar abundance. The resulting lEM profile, along with the 
best fitting model, is shown in Figure \13\ The deviant 
point at 2 kpc is due to small-scale clunipiness in the 



ICM, visible in the X-ray image. We find a density jump 
of 1.97 ± 0.12, which corresponds to a M = 1.71 ± 0.1 
shock with a temperature jump of 1.72 ± 0.12. Unfortu- 
nately, there are too few counts to derive a deprojected 
temperature profile to compare to the inferred temper- 
ature jump, although the corresponding inner jump in 
the azimut hall y averaged temperature profile of about 
1.1 (Figure 10 1 is much smaller than the jump of 1.7 pre- 
dicted from the density discontinuity. To better quantify 
the temperature jump across the 1.5 kpc shock, we ex- 
tracted a spectrum of the total emission from the post- 
shock region in the core wedge used to fit the surface 
brightness profile, and from a similar sized region in the 
same wedge just outside the shock edge. We find that 
the (inner) post-shock region is best fit by a two temper- 
ature thermal model, with best-fitting temperatures of 
kTi = 0.641°:°;^ and fcTa = 0.97j:!]:^^, while the pre-shock 
region is adequately described by a single thermal model 
with kT = 0.5 ±0.01. If we assume that the two thermal 
components in the post-shock region give the tempera- 
tures of the shocked gas and local projected pre-shocked 
gas this gives a temperature jump of 1.5 ± 0.3. This is 
consistent with a temperature jump of 1.7 as expected 
from the Mach number, although the large uncertainties 
limit the usefulness of this comparison. The properties of 
the 10 kpc and 1.5 kpc shocks are summarized in Table[2| 
For both the 1.5 kpc and 10 kpc shocks, the observed 
temperature rise is less than what is expected based on 
the Mach number, even in the deprojected temperature 
profiles. To better estimate the expected measured tem- 
perature rise associated with the shocks we ran ID hydro- 
dynamical simulations to fully model their evolution. For 
the 10 kpc shocks, the simulations start with an isother- 
mal sphere with kT — 0.67 and a power law density pro- 
file, in hydrostatic equilibrium. The logarithmic slope of 
the density profile was taken to be 1.46, which was deter- 
mined by fitting the slope of the surface brightness pro- 
file beyond the 10 kpc shocks. The shocks were initiated 
by a central point explosion and allowed to propagate 
freely. Since the gas is initially isothermal with a power 
law density profile, the simulations are scale free, so we 
can choose the point of evolution at which the code best 
reproduces the observed surface brightness jump at the 
shock front and scale accordingly. The resulting Mach 
number for the 10 kpc shocks is M — 1.52, in excellent 
agreement with what was found above by fitting a discon- 
tinuous power law density model to the integrated emis- 
sion measure profile. The expected projected tempera- 
ture profile is shown in Figure |14[ The solid lines show 
the emission weighted temperature (with 90% confidence 
ranges) and the points show single temperature fits to 
more accurately simulated model spectra folded through 
the Chandra response (the discrepancy between emission 
weighted temperatures and fitted temperatures for multi- 
temperature gas has previously been noted by Mazzotta 
et al. 2004). The simulations predict a projected temper- 
ature rise of ^ 0.1 keV, consistent with observations (see 
Figure [lOl). Similar simulations fit to the 1.5 kpc shock to 
the souEneast predict a measured projected temperature 
rise of ~ 30%, consistent with observations (although the 
observational uncertainties are large). We conclude that 
the projected temperature profiles are consistent with 
the calculated Mach numbers (although the temperature 
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structure of the 1.5 kpc shock is poorly constrained). 
The deprojected temperature measurements cannot re- 
solve the temperature jump due to the narrow width of 
the shock and the rarefied cool gas behind the shock (the 
temperature of which drops slightly below the ambient 
temperature in our isothermal simulations). 

5.2. The Outer Edge at 25 kpc 

There is a weak surface brightness edge (most clearly 
seen in Figure|3| surrounding most of NGC 5813 at a dis- 
tance of ^160 (25 kpc) from the central AGN. To check 
the significance of this feature we extracted the surface 
brightness profile across this edge in two sectors, one to 
the northwest between 12° - 79° (measured north from 
west, the same angular range spanned by the bright sec- 
tion of the elliptical edge to the northwest) and a wider 
wedge between 206° - 320° to the south. The wedges 
were chosen to match the curvature of the outer edge, 
so that the profiles were off center from the peak of the 
overall diffuse emission. Th e re sulting surface brightness 
profiles are shown in Figure 



15 



Both profiles show a sim- 
ilar change in slope at ^ 170^^ (26 kpc), at the position 
of the outer edge, and the northwestern profile shows a 
sharp jump at the same location. We conclude that the 
outer edge-like feature indicated in Figure [3] is real, and 
corresponds to a change in slope and possibly a discon- 
tinuity in the surface brightness profile. The edge lies 
just beyond the outer X-ray cavities indicated in Fig- 
ure [3] and the association is reminiscent of the 1.5 kpc 
and 10 kpc shocks to the inner and intermediate cavities. 
The discontinuous jump is also stronger in the northwest 
than in the south, consistent with what is seen for the 
elliptical 10 kpc shock front edge. It may therefore rep- 
resent the weak remnants of a shock associated with this 
older outburst. 

To determine the nature of this edge, we fit the surface 
brightness profile in the wedge to the the south, where 
this feature is the sharpest. Although there is a signif- 
icant surface brightness discontinuity to the northwest 
(see Figure 15 1, we focus on the wedge to the south, since 
we have better statistics in the wider southern wedge, 
and since the shape of the northwestern discontinuity 
is not well-described by our shock model density pro- 
file (compare the profile shapes in Figure 15 and in Fig- 



ure [12). The shape of the surface brightness profile is 
obviously not well-described by a single power-law, and 
modeling this profile with a projected power-law density 
model did not provide an acceptable fit. We also fit the 
profile with a /3-model density profile, which gave a poor 
fit with = 13.9/6. The fitted model showed asymmet- 
ric residuals in the region of the surface brightness edge, 
so we tried fitting the edge with a proj ecte d discontin- 
uous power-law density model, as in § |5.1| This gave 



an improved fit, with xt — 0-4/4, and a density jump 
factor of ^outer = 1.06^Q Jg, consistent with no jump. 
The inner and outer slopes of the density profile were 



-0.62 



-1-0.09 
-0.18 



and P = -2.30 



+0.25 
-0.32- 



Unfortunately, 



there were inadequate statistics to measure the temper- 
ature and abundance profiles across this outer edge and 
confirm it as an old shock. In particular, the edge could 
in principle be a metallicity edge, or the interface be- 
tween the galaxy atmosphere and the extended group 
atmosphere (where one would expect to see a change in 



slope of the surface brightness profile). We conclude that, 
while the outer edge is consistent with an old shock with 
M « 1.1 associated with the outermost X-ray cavities, 
the data are also consistent with no density jump and 
further observations are needed to determine its nature. 

5.3. The X-ray Cavities 

The X-ray image shows three pairs of roughly collinear 
cavities (Figure fl| Figure |4]) . Although most of the cav- 
ities are clear, the outermost cavities (in particular, the 
southwestern outer cavity) are weak features. To check 
their significance, we extracted surface brightness profiles 
across each outer cavity, centered on the overall di ffuse 
emission. The resulting profiles are shown in Figure 16 
Each profile shows a significant dip at the location of me 
outer cavities. The profile across the northeastern cav- 
ity also shows a significant hump just beyond the cavity 
dip, corresponding to the bright outer rim seen in the X- 
ray images. For the southwestern outer cavity, the rise 
just outside the cavity is less pronounced, since it lacks a 
bright rim. We conclude that the faint outer southwest- 
ern cavity indicated in Figure [3] is likely a real feature 
and represents a paired cavity to the outer northeastern 
cavity, each initially created by the same AGN outburst 
from the central SMBH. 

Each of the three pairs of X-ray cavities in NGC 5813 
is likely associated with a distinct AGN outburst. Al- 
though the position angle of the inner and intermediate 
cavity pairs appears to vary slightly, by about 10° - 15°, 
possibly indicating that the central black hole that has 
inflated the cavities is precessing, the cavities are roughly 
collinear. The fact that the cavities are regular and 
collinear suggest that they have evolved passively, i.e., 
have not been disturbed by gas motions due to sloshing, 
turbulence, mergers, etc., consistent with previous con- 
clusions on the dynamical state of this system (Tran et 
al. 2001; Emsellem et al. 2007). Measuring the proper- 
ties of these cavities is therefore of great interest, since 
a comparison of the different pairs will provide informa- 
tion on how the cavities evolve, and on the outburst his- 
tory of the central AGN (in particular, we want to know 
whether the outbursts have a similar total energy output 
and mean power, as might be expected in a near "steady 
state" AGN feedback model, or whether the outbursts 
change signiflcantly even in an otherwise apparently re- 
laxed system hke NGC 5813). 

5.3.1. Cavity Ages and Energies 

Table [3] summarizes the properties of the X-ray cav- 
ities, which are shown as overlaid regions in Figure [4j 
We assume that the cavities rise in the plane of the sky, 
that the mass of material within the cavities is negligible, 
and that each cavity has a spherical or oblate spheroidal 
geometry with the minor axis in the plane of the sky. 
The former assumption is suggested by the fact that the 
cavities are detected in the image, and that the inner 
two pairs are near the same projected radii as their as- 
sociated shock fronts. If the global gas distribution was 
significantly extended along the line of sight, with the 
bubbles rising along this extension, they would be diffi- 
cult to detect for large inclination angles of the extension 
due to the large column of cool gas at smaller radii be- 
hind (or in front of) the bubbles (see Birzan et al. 2009). 
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Columns 2 & 3 in Table [3] give the cavity major and mi- 
nor axis, respectively, ana Column 4 gives the distance 
from the AGN to the center of the cavity. The rise time 
of the cavities, given in Column 5, is calculated assum- 
ing that they rise buoyantly at half the sound speed , 
similar to what is found from simulations (e.g., 0.6-0.7cs 
in Churazov et al. 2001). For a kT = 0.65 keV gas, the 
sound speed is Cg — 416 km s~^. Note that for the inner- 
most cavities, the distance from the AGN is comparable 
to the cavity size, suggesting that they are currently be- 
ing, or have only recently been, inflated by the AGN. 
Therefore the rise times, which are computed using the 
distance from the AGN to the cavity center, are not a 
reliable estimate of the true age, since early on the cav- 
ities are driven by the momentum of the jet (see Birzan 
et al. 2008). 

Finally, the mechanical energy required to inflate the 
bubbles is given in Column 6, which we estimate as PV , 
where P is the pressure at the location of the bubble cen- 
ter (taken from the azimuthally averaged pressure profile 
shown in Figure 10 ) and V is the cavity volume. The to- 



tal internal energy of each cavity is expected to be a few 
times the mechanical energy (^ 3PV for a relativistic 
plasma, see McNamara & Nulsen 2007). 

The total mechanical energy is about the same for the 
southwestern and northeastern cavities for both the in- 
nermost (~ 1.3 X 10^^ erg) and intermediate (~ 1.4 x 
Iq56 gj-g-j pa^jj-g (taking the sum of the two intermedi- 
ate cavities to the northeast), consistent with each pair 
having formed in two distinct AGN outbursts. The 
mechanical energies for the outer cavities differ signifi- 
cantly from one another (2.6 x 10^^ erg versus 6.0 x 10^^ 
for the northeastern and southwestern cavities, respec- 
tively). We suggest three possibilities to explain this dis- 
crepancy. First, the outer southwestern cavity is only 
marginally detected, so its measured properties may be 
inaccurate. Second, the outer cavities may be in the pro- 
cess of breaking apart, as suggested by their significantly 
different volumes, and if they are not devoid of X-ray 
emitting gas, then the mechanical energy will be incor- 
rectly estimated. Finally, Figure [l] shows that the inter- 
mediate northeastern cavities may connect to the outer 
northeastern cavity. If this is indeed the case, energy may 
"leak" from the intermediate to the outer cavity, adding 
to its internal energy and inflating it, making it easier to 
detect (there is no such connection in the southwest). 

A comparison of the ages of the cavities (approximated 
by the rise times given in Table J3| and the shocks (ap- 
proximated by the travel time given in Table [2| for the 
"inner" and "middle" features shows that they are sim- 
ilar (^ 10^ yr), consistent with the interpretation that 
each set of features was formed at the same time by the 
same outburst event. The cavity ages are systematically 
larger than the shock ages derived from our hydrodynam- 
ical simulations (by about a factor of 3). This is likely 
due to the cavities being initially driven and inflated at 
some significant distance from the central AGN, only to 
rise buoyantly after the end of the outburst, so that re- 
garding the cavities as buoyantly rising bubbles early in 
their lives is not accurate. Thus, using the buoyant rise 
velocities of the cavities likely overestimates their ages. 

5.3.2. Pressure Balance with the ICM 



The non-thermal radio pressure in X-ray cavities, un- 
der the assumptions of hydrostatic equilibrium and elec- 
tron dominated pressure, is commonly found to be less 
than the pressure of the surrounding gas derived from 
X-ray observations (Birzan et al. 2008). To accurately 
estimate the radio pressure, flux measurements at mul- 
tiple frequencies are required to estimate the spectral 
index ar of the radio emission. In NGC 5813, only the 
innermost cavities are detected at more than one fre- 
quency (see § 3.1). Unfortunately, the large beam size 



and contamination from extended emission outside the 
cavities at 235 MHz make it difficult to accurately mea- 
sure ar- Comparing the total emission in the region of 
the inner cavities, including contributions from extended 
emission outside of the central cavities at 235 MHz and 
1.36 GHz (the latter detected in archival VLA C array 
configuration observations) and emission from the core, 
gives — 0.88. If we assume — 0.7 within the inner 
cavities, a typical value for young radio lobes (Birzan et 
al. 2008), we find that the non-thermal pressure in each 
cavity is Prad ^ 2.3 x 10^^^ erg cm~'^ (where we have 
used the revised equipartition equations of Brunetti et 
al. 1997 with a low-energy cut-off of 7min = 100). This 
is much less than the X-ray gas pressure outside the cav- 



ities P, 



gas 



10- 



-10 



erg cm (see Figure 10 1. Balancing 



the pressure with low energy electrons alone would re- 
quire Qfj, sa 2, which is unusually steep for young radio 
cavities (Birzan et al. 2008). If we assume = 0.7 and 
vary the ratio of the energy in protons to the energy in 
electrons k, we find that k w 2000 is required to balance 
the non-thermal and X-ray gas pressures, which is typ- 
ical for radio galaxies in cool cores (Birzan et al. 2008; 
Dunn et al. 2010; Gitti et al. 2010). 

We compared the radio and cavity power of the in- 
ner cavities with predictions based on the Birzan et al. 
(2008) sample. The 1.36 GHz radio power for each of 
the inner cavities is Piaeo ~ 1 x 10^° W Hz~^. The re- 
lation of Birzan et al. (2008) then predicts a total cavity 



power of 



9 X 10 erg s . Taking the cavity 



ages and mechanical energy (PV) from Table [3] and as- 
suming total cavity energies of 3PV gives a power of 
Pcav ~ 4 X 10'''^ erg s~^ for the inner cavities, more than 
an order of magnitude less than the estimate from Birzan 
et al. (2008). However, it is within the range of the large 
scatter in the Birzan et al. (2008) sample (see their Fig- 
ure 6). More recently, Cavagnolo et al. (2010) derive a 
relation between Pcav and Pradio for lower mass systems 
and find a steeper slope than Birzan et al. (2008) . Their 
sample includes results from the Chandra observations of 
NGC 5813 we consider here, so that NGC 5813 is con- 
sistent with their derived relation. 

5.4. Buoyantly Lifted Gas 

The temperature map (Figure p| shows an extension 
of cool gas along the line defined oy the X-ray cavities, 
offset ^ 3 kpc to the southeast. This cool gas extends 
11 kpc to the inner edges of the intermediate cavities. 
The most natural interpretation of this feature is cool 
gas that has been buoyantly lifted by the intermediate 
X-ray cavities, as seen in other systems (Forman et al. 
2007; Simionescu et al. 2008; Kirkpatrick et al. 2009; 
Simionescu et al. 2009; Randall et al. 2010; Werner et al. 
2010). The Ha image, taken with the SOAR telescope 
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(Donahue et al. 2007), with the X-ray surface bright- 
ness contours overlaid, as well as the temperature map 



with the Ha contours overlaid, are shown in Figure 17 



The Ha filaments are co-spatial with the trail of cool 
gas seen in the X-ray temperature map, as seen in many 
other systems (e.g., Sanders et al. 2007, 2009), confirm- 
ing the presence of cool gas. Assuming line ratios of 
N[n]A6583/Ha = 1 and N[H]A6548/Ha = 0.35, the total 
Ha flux is Fiia = 9.2 x lO"^'* erg s^^ cm^^ (assuming 
no intrinsic absorption). The Kennicutt relation for the 
star formation rate 



SFR(M0/yr) = 7.9 x W'^^Lnaierg s" 



(1) 



(Kennicutt 1998) then gives SFR = O.O9Af0 yr^^. This 
assumes that the Ha emission is completely driven by 
UV radiation from young stars, which may not be the 
case (e.g., Fcrland et al. 2009), and that there is no in- 
trinsic absorption. A violation of the former assump- 
tion would give an overestimate of the SFR, while a vi- 
olation of the latter would give an underestimate. We 
compared the derived SFR to the mass cooling rate es- 
timated from the X-ray observations by fitting the spec- 
trum of the cool filament with a vapec plus a cooling 
flow VMCFLOW model, with the abundances fixed at the 
best fittin g values for the total diffuse emission given 
in § |4.2.4 This gives a mass cooling rate in the gas 
of 0.41 ± 0.07 yr^^. One concern with this result is 
that fitting a multi-temperature gas with a single cooling 
flow component may boost the inferred mass accretion 
rate. Unfortunately, we were unable to obtain a fit with 
reasonable parameter constraints from the Chandra data 
using a model combining an APEC and two VMCFLOW 
components. We therefore consider the mass accretion 
rate inferred from the fit to XMM-Newton RGS data de- 
scribed in 



4.2.4 



From the two cooling fiow compo- 
nents, we infer upper limits on the mass accretion rate 
of < 0.45 Mq yr^^ and < 0.25 M© yr~^ for gas cool- 
ing above and below 0.4 keV, respectively (these are 2cr 
upper limits). Therefore, we conclude that the upper 
limit on the mass accretion rate from X-ray observations 
(< 0.25 Mq yr~^) and the star formation rate implied 
by Ha observations (SFR = O.O9A/0 yr~^) are consistent 
with star formation being fueled by gas cooling down 
from X-ray temperatures. 

For the cool gas to be buoyantly lifted by the X-ray 
cavities, its total mass must be less than the mass of 
gas displaced by the cavities. In particular, simulations 
indicate that the mass of gas buoyantly lifted by an 
AGN-blown bubble is about half the mass displaced by 
the bubble (Pope et al. 2010). We estimated the gas 
mass in the southern filament by fitting the spectra in 
this region with an APEC model and assuming that the 
cool gas is contained in a cylinder of radius 3.5 kpc and 
length 9.2 kpc, with the axis in the plane of the sky. 
The emission is dominated by the cool gas in this region, 
and accounting for the projected hot gas did not signif- 
icantly change our results. We find a total gas mass of 
Mfii = 1.5 X 10* M0. A similar fit to an annulus sur- 
rounding the southwestern intermediate cavity gives an 
average electron density of Ue = 0.022 cm~'^, giving a 
total mass of displaced gas Mdisp = 1-6 x IO^A/q, similar 
to the mass of gas in the filament. Thus, Mju is larger 
than the value predicted by simulations 0.5Mdisp) 



by a factor of two. We note that deep observations of 
the buoyantly lifted filaments in M87 reveal that they 
have a fine filamentary structure (Forman et al. 2007), 
in contrast with the solid cylindrical geometry we have 
assumed above (Werner et al. 2010 argue that the filling 
factor in M87 is of order unity in most regions, although 
they find a filling factor less than unity in some regions 
where there is fine filamentary structure). If the filling 
factor is less than unity then Mfu will be smaller by the 
same fraction. We conclude that the gas mass of the 
cool filament to the south is consistent with having been 
buoyantly lifted by the intermediate southwestern cavity, 
and is consistent with simulations if the filling factor is on 
the order of ^ 0.5. We also note that if the filament has 
indeed been buoyantly lifted by the intermediate cavity 
then the filament, and hence the trajectory of the inter- 
mediate cavity, cannot lie far from the plane of the sky 
without assuming a small filling factor, since the length 
of the filament (and hence the volume it occupies) grows 
with inclination angle. 

5.5. The Offset of the Central AGN 

As noted in § [3] the central AGN is offset ~ 0.5 kpc 
southeast of the line defined by the X-ray cavities (see 
Figure l2|. It is also offset from the center points of the 
elliptical edges defined by both the 1.5 kpc and 10 kpc 
shock fronts (which are each roughly coincident with the 
line defined by their respective cavity pairs), by about 
2.3" (400 pc) for the 1.5 kpc shock and 7.5" (1.15 kpc) 
for the 10 kpc shock. A comparison of the Chandra 
image with the optical Sloan Digital Sky Survey image 
(Adelman-McCarthy et al. 2008) shows that the centroid 
of the optical emission is coincident with the AGN, and 
is separated from the center points of the elliptical X- 
ray shock fronts. This suggests that the cD galaxy has 
some peculiar velocity relative to the ICM, and that the 
AGN has moved since first inflating the X-ray cavities. 
The initial outburst that inflated the inner cavities oc- 
curred about 3 X 10^ yr ago (see Table [2]). Since then, 
the AGN has traveled about 400 pc in projection. This 
requires a relative velocity between the galaxy and the 
ICM of at least ~ 130 km s~^(this is a lower limit since 
we measure the projected velocity). Similarly, the 10 kpc 
shock implies a relative velocity of ^ 100 km s~^. This is 
smaller than the host group's (NGC 5846) velocity dis- 
persion (322 km s-\ Mahdavi et al. 2005), and the rel- 
ative radial velocity between NGC 5813 and NGC 5846 
260 km s~^), and is therefore a reasonable peculiar 
velocity for the cD galaxy relative to the group mean 
(or for the flow velocity of the ICM gas). We suggest 
that during the last ^ 2 x 10^ yr there has been a rela- 
tive motion between the NGC 5813 galaxy and its ICM, 
from northwest to southeast, in projection, with the cen- 
tral AGN driving outbursts (i.e., inflating cavities into 
the ICM) at two locations (the outer cavities and edge 
are too poorly resolved to accurately measure the loca- 
tion of the associated outburst). The relative velocity 
may be due to the peculiar velocity of the cD relative to 
the group mean, or to bulk gas motions or gas sloshing 
of the ICM, or both. We note that once the outburst 
begins and the cavities are initially inflated, the location 
at which the energy is injected into the cavities is unim- 
portant. The AGN may continue to move off of the axis 
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of symmetry of a pair of cavities while still depositing en- 
ergy into them, driving their expansion and the resulting 
shocks. 

5.6. Outburst Energy 

During an AGN outburst, the jets inflate cavities in 
the ICM, which do work on the surrounding gas and 
drive shocks. It is convenient to regard the work done 
by the expanding lobes as "shock energy" , in which case 
the outburst energy that is available to heat the ICM 
is deposited to the ICM in two forms: the internal en- 
ergy of the X-ray cavities and the shock energy. Since 
NGC 5813 shows both cavities and shocks from two dis- 
tinct outbursts, we can compare the total energy, mean 
power, and the energy budget between shocks and cavi- 
ties for each outburst. The shock energy can be roughly 
estimated from the pressure increase across the shock 
front. If a total energy E is added to a gas of volume 
V the pressure increase is roughly AP ^ E/V. For a 
shock with a known Mach number, the ratio of the post- 
and pre-shock pressure fp = {P + AP) /P is given by 
the Rankine-Hugoniot shock jump conditions. The total 
shock energy is therefore 



E « PVifp - 1). 



(2) 



As a consistency check, we compare the estimated shock 
age, total energy, and mean power with results from our 
hydrodynamical model in Table [2j The shock age iagc.cst 
is estimated as the travel time from the current position 
of the shock front to the center point of the elliptical 
shock edge, assuming a constant Mach number and a 
sound speed of Cs = 416 km s~^ for a, kT = 0.65 keV 
gas. The shock energy iJsh.cst is estimated using equa- 
tion [2] assuming a prolate ellipsoidal geometry for the 
volume contained within the shock front with semi-major 
and -minor axes of 1.44' (13.5 kpc) and 1.13' (10.5 kpc) 
for the 10 kpc shocks and 17.6" (2.7 kpc) and 11.7" 
(1.8 kpc) for the 1.5 kpc shocks. Pressures were taken 
from the aziniuthal pressure profile shown in Figure |10[ 
Although the energies we give from the hydrodynami- 
cal model represent the total outburst energy, the point 
explosion model minimizes the internal energy in cavi- 
ties, such that the energy in the central cavity is only a 
few percent of the total outburst energy. The model en- 
ergy is therefore a very good approximation of the shock 
energy, and we refer to it as such. The model shock en- 
ergies i?sh, model were scaled to account for the difference 
in total volume for the spherically symmetric model as 
compared to the observed elliptical edges, and for the 
lower average pressure along the elliptical shock fronts, 
which cover a range in radii, assuming the shock energy 
scales as in equation [2] (these effects somewhat balance 
each other, as the volume correction increases the total 
energy, while the pressure correction decreases it). The 
correction factors are between 15-40%. The estimated 
and model shock energies and mean powers agree rea- 
sonably well, within factors of a few, demonstrating the 
consistency between rough estimates and results from our 
point explosion hydrodynamical model. 

Table [2] indicates that the shock energy for the cur- 
rent outburst is more than an order of magnitude smaller 
than for the previous outburst. The mean power of the 
current outburst is also less than that of the previous 



outburst, by about a factor of six (1.5x10''^ erg s ^ ver- 
sus 1.0x10''^ erg s~^, where we take the mean outburst 
power to be the sum of the shock energy and the 3PV in- 
ternal energy of the cavities divided by the shock model 
ages). The lower shock energy of the most recent out- 
burst may indicate that it is ongoing, having only de- 
posited ^1/40 of its expected total energy output into 
the observed shocks (assuming that the current outburst 
is similar to the previous one). However, the lower shock 
energy may simply be a result of the lower mean power 
of the current outburst. We note that from X-ray ob- 
servations of elliptical galaxies Allen et al. (2006) find 
evidence that accretion flows around central AGN are 
stable over a few million years, whereas we find that the 
mean jet power varies on time scales of ~ 10^ yr (the 
time between outbursts). 

We wish to compare the energy in shocks to the en- 
ergy in the X-ray cavities for each outburst. The to- 
tal internal energy of the cavities is roughly 3 times the 
mechanical energy (PV) required to inflate the cavities 
(McNamara & Nulsen 2007), which is given in Table [3| 
We find total cavity internal energies of 8.6 x 10^^ erg 
and 7.8 x lO^'^ erg for the outbursts that produced the 
intermediate and inner cavities, respectively. Thus, the 
total internal energy in cavities is roughly 30% of the 
shock energy for the previous outburst and 1.3 times the 
shock energy for the current outburst. This is consistent 
with the current outburst being young, with the jets ac- 
tively inflating the inner cavities and driving the 1.5 kpc 
shock, whereas for the previous outburst the shock has 
detached from the cavities, which rise buoyantly and lose 
energy as they age. We note that while our point explo- 
sion shock model is only approximate, so that the shock 
energies are somewhat uncertain, the relative sense of the 
energies is correct. Therefore, the larger fraction of total 
energy in cavities in the current outburst (as compared 
to the previous one) is a robust result. 

For the outermost pair of cavities, the total internal 
energy in the cavities is ^ 9.6 x 10^^ erg, on the order of 
the total energy of the intermediate cavities. However, 
80% of this energy is contained in the northeastern cavity. 
While it is in principle possible for the mechanical energy 
to differ between paired cavities from the same outburst, 
the measurements for the southwestern cavity are uncer- 
tain since it is only marginally detected. Furthermore, 
such a large difference is not seen in the inn er and in- 
termediate cavity pairs. As discussed in § |5.3[ the outer 
cavities may be in the process of breaking apart, and the 
outer northeastern cavity may have been re-energized by 
the intermediate cavity, making the measured mechani- 
cal energy (and hence the internal energy) uncertain. 

We conclude that the lower total energy of the most 
recent outburst, as compared to the previous outburst, 
suggests that the outburst is ongoing, though this may 
simply be an effect of the lower mean power. Al though 



the luminosity of the central source is low (see § 4.2.5) 
it may be in a short term quiescent state or heavily 
obscured, and does not rule out an ongoing outburst. 
AGN luminosities are known to vary by several orders 
of magnitude over very short time scales (e.g., Harris et 
al. 2009). Results from our ID hydrodynamical simula- 
tions suggest that the mean power over longer timescales 
(^ 10^ yr) can also vary signiflcantly between outbursts. 
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even in an otherwise relaxed system like NGC 5813. 

5.7. The Balance Between Heating and Cooling in the 

ICM 



As noted in § 5.6 the total mechanical energy output 
of the AGN that is available to heat the ICM is primar- 
ily in two forms: the internal energy of X-ray cavities, 
which rise buoyantly after being inflated by jets from the 
AGN, and the "shock energy" , from shocks driven by the 
rapid inflation of the cavities. When and where the in- 
ternal energy of the X-ray cavities gets transferred to the 
ICM is not well understood. In contrast, the local heat 
input at the shock front can be calculated directly from 
the Mach number. Furthermore, shock heating has two 
desirable features. First, the gas is more strongly heated 
in the core where Mach numbers are larger, close to the 
central AGN, which is the region of interest for regulat- 
ing feedback between the ICM and the central SMBH. 
Second, the heating is roughly isotropic (as opposed to 
heating with jets or with the internal energy of the X-ray 
cavities). We therefore consid er th e energy input due to 
shocks, which, as we show in § |5.6[ contain 40-80% of the 
total outburst energy in NGC 5813 (for a detailed dis- 
cussion of heating with AGN outburst shocks see David 
et al. 2001). 

To balance radiative cooling with AGN feedback 
shocks, the average outburst power must be on the order 
of or larger than the rate of radiative cooling. We esti- 
mated the radiative cooling rate as the X-ray luminosity, 
which we obtained by fitting the spectrum from the total 
emission within 170" (26.3 kpc) with an absorbed APEC 
model (using a VAPEC or two APEC model did not sig- 
nificantly change the resulting luminosity). The derived 
X-ray luminosity of Lx = 5.4 x lO''^ 
implies a mean cooling time of t 
take the cooling time to be the time it would take to 
radiate away all of the internal energy of the gas at the 
current luminosity). For the 10 kpc shocks, we calculate 
the total shock energy to be ^ 3 x 10^'' erg (see § 5.6). 
Therefore, only 6 such outbursts are needed per cooling 
time to provide enough energy to completely offset ra- 
diative cooling with shocks alone within 26 kpc. The 
time between outbursts given by both the buoyant rise 
times of the cavities (Table |3 1 and, more reliably, from 
hydrodynamical simulations of the shocks (Table pi) is on 
the order of 10^ yr, allowing 100 shocks per cooling time, 
more than is needed to provide the necessary energy. 

While results from the 10 kpc shocks indicate that 
there is sufficient energy in shocks alone to offset radia- 
tive cooling, there are two points that must be consid- 
ered. First, as discussed in § 5.6 the observed shocks 
suggest that the current and previous outburst shock en- 
ergies differ in strength by more than an order of mag- 
nitude (although the most recent outburst may be ongo- 
ing). Furthermore, only some fraction of the total shock 
energy will go into heating the gas within the cooling ra- 
dius, and for weak shocks this fraction is relatively small 
(< 10%). We therefore consider the /omZ balance of shock 
heating and radiative cooling for each shock. The pri- 
mary effect of radiative cooling is to reduce the entropy 
of the gas. To offset cooling, the heating mechanism is 
required to increase the entropy by at least this amount. 
Shocks will also affect the kinetic, thermal, and potential 



erg s ^ within 26 kpc 
cool = 1.0 Gyr (where we 



energy of the gas, but these effects are transient, and for 
heating it is the change in entropy of the gas that is rele- 
vant. Hence, the entropy increase AS* caused by a weak 
shock can offset a radiative heat loss of AQ — T AS, 
where T is the gas temperature. Expressed as a fraction 
of the gas thermal energy, the effective heat input from 
one shock is therefore {TAS)/E = Aln(P/p^), where 
E is the thermal energy and 7 is the adiabatic index. 
For the 10 kpc shocks, the Mach numbers given in Ta- 
blel2]imply a change in \n{P/p'^) across the shock front of 
^ 5%. Therefore, to balance the total entropy decrease 
of the gas about 1/0.05 = 20 outbursts are needed per 
local cooling time to completely offset cooling with shock 
heating. By the same argument, 10 outbursts per cooling 
time are needed to replace the thermal energy of the gas 
at the 1.5 kpc shock. The cooling time of the pre-shock 
gas just outside the 10 kpc shock edge is 9 x 10^ yr, so 
an outburst interval of 10^ yr gives 90 shocks per local 
cooling time. For the 1.5 kpc shock, the pre-shock gas 
cooling time is 2 x 10® yr, giving 20 shocks per cool- 
ing time. Thus, we conclude that shock heating alone is 
sufficient to offset radiative cooling of the gas within the 
1.5 kpc and 10 kpc shock fronts, consistent with previous 
suggestions that shocks can offset cooling close to central 
AGN (Nulsen et al. 2007). The X-ray cavities rise buoy- 
antly, and release their internal energy to heat the ICM 
gas at larger radii. Although previous studies have found 
other systems where there is enough total shock energy 
to offset radiative cooling (e.g., M87 Forman et al. 2005, 
with M87 also showing concentric shock fronts from mul- 
tiple AGN outbursts; Hydra A Nulsen et al. 2005), here 
we explicitly show that the fraction of shock energy that 
goes into heating the gas (5-10%) is sufficient to balance 
radiative cooling locally at the shock fronts. The out- 
burst interval is short enough for such shocks to offset 
cooling over much longer timescales. This demonstrates 
that AGN feedback can operate to heat the gas within 
galaxies, as well as the more extended ICM in clusters, 
as required by some galaxy formation models (e.g., Kor- 
mendy et al. 2009). 

6. SUMMARY 

We have presented results based on Chandra, VLA, 
GMRT, and SOAR observations of NGC 5813, the dom- 
inant member of a galaxy group. The ICM shows 
clear signatures from three distinct AGN outbursts 
in an otherwise relaxed system (including two shocks 
with detectable temperature jumps), making this object 
uniquely well-suited to the study of AGN feedback. We 
find the following: 

1. Three pairs of collinear cavities, where each pair 
is associated with a distinct AGN outburst. The 
inner two pairs are associated with unambiguous 
shocks (with Mach numbers Mi = 1.7, Mo = 1.5), 
with clear temperature rises, that can be directly 
detected from the X-ray data. The outermost cav- 
ity pair also has an associated surface brightness 
edge. The properties of this edge are consistent 
with an old shock associated with the outermost 
cavities, although the current data do not rule out 
other interpretations (e.g., a transition region from 
the galactic atmosphere to the extended group at- 
mosphere). The locations of the cavities and the 
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shocks indicate an outburst interval of ^ 10^ yrs. 

2. Diffuse radio emission, filling the inner cavities at 
1.36 GHz and 235 MHz and the intermediate cav- 
ities at 235 MHz. Radio emission from the outer 
cavities is not detected. This reflects the greater 
age of the relativistic particles in more distant cav- 
ities. 

3. A cool trail of gas that has been buoyantly lifted 
by the intermediate cavities. The gas is co-spatial 
with filaments seen in Ha observations. 

4. The mean power of the current outburst is six times 
less than that of the previous outburst, indicating 
that the average jet power can vary significantly be- 
tween outbursts, even in near "steady state" AGN 
feedback with regular outbursts in an otherwise dy- 
namically relaxed system. 

5. The total heat energy input from shocks alone is 
sufficient to balance radiative cooling locally. This 
heating takes place in the core, close to the central 
AGN, which is the region of interest for regulating 



feedback between the ICM and the central SMBH. 
Thus, for the first time, we explicitly show a sys- 
tem where shock heating alone can locally balance 
radiative cooling and regulate AGN feedback. 
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TABLE 1 

Summary of the Radio Observations 



Radio 


Project 


Observation 


Array 


Frequency 


Bandwidth 


Integration 


HPBW, PA 


rms 


telescope 




date 




(MHz) 


(MHz) 


time (min) 


("x", °) 


(/.Jy b-l) 


GMRT 


14SGA01 


Aug 2008 


full 


235 


8 


100 


16.5 X 15.1,74 


300 


VLA 


AF0188 


Apr 1990 


A 


1490 


50 


45 


1.3 X 1.2,28 


22 


VLA 


AW0202 


Jan 1988 


B 


1360 


50 


30 


4.9 X 4.7,-84 


20 


VLA 


AC0488 


Sept 1997 


C 


1360 


50 


6 


27.4 X 16.8,53 


30 


VLA 


AW0112 


Jun 1984 


C 


4860 


50 


8 


5.0 X 4.5, -4 


25 



TABLE 2 
Properties of the Shocks 



ID 


(kpc) 






X d 
''age, model 

(10^ yr) 


-'-'shsmodel 

(1057 erg) 


l^model' 

(10^2 erg s-1) 


MaeeS 
(IO^Mq) 


/ h 
'-age, est 

(107 yr) 


-^sh,cst 

(1057 erg) 


(lO*'^ erg s-i) 


Inner, SE 


1.4 


-0.12 


1.71 


0.3 


0.06 


0.6 


0.3 


0.2 


0.2 


2.7 


Middle, SE 


9.9 


1 

J^-O^-0.07 


1.48 


1.4 


2.2 


5.0 


8.2 


1.6 


3.0 


5.9 


Middle, NW 


11.5 


, 7c:+0.04 


1.53 


1.3 


3.3 


8.0 


22.4 


1.6 


3.0 


5.9 



^ Distance from the AGN to the shock front. 
Density jump at shock front. 
Mach number. 

Shock age, from the hydrodynamical model. The point explosion approximation gives large Mach numbers at early times, so the shock 
ages are underestimated. 

° Shock energy, from the hydrodynamical model. 

' Mean shock power, from the hydrodynamical model. 

s Minimum accreted mass needed to power the model outburst (assuming 100% efficiency). 

Shock age, estimated as the travel time from the current position to the center point of the elliptical shock edge. Since the current Mach 
number is used to give a constant velocity, the shock ages are overestimated. 
' Shock energy, estimated using Equation l2| 
J Estimated mean shock power. 



TABLE 3 
Properties of the X-ray Cavities 



ID 




fob 






^mech 




(kpc) 


(kpc) 


(kpc) 


(107 yr) 


(1055 erg) 


Inner, SW 


0.95 


0.95 


1.3 


0.6' 


1.1 


Inner, NE 


1.03 


0.93 


1.4 


0.7' 


1.5 


Middle, SW 


3.9 


3.9 


7.7 


3.6 


15.3 


Middle-lf5, NE 


2.9 


2.2 


4.9 


2.3 


9.3 


Middle-2f5, NE 


2.8 


2.4 


9.3 


4.4 


4.1 


Outerh, SW 


5.2 


3.0 


22.2 


10.4 


6.0 


Outer, NE 


8.0 


4.4 


18.0 


8.5 


26.0 



Semi-major axis. 
^ Semi-minor axis. 
Distance from central AGN. 

Lower limit on the bubble rise time, assuming that each bubble rises at half the sound speed for a 0.65 keV gas. 
° PV mechanical energy required to inflate the cavity. 

' The cavity size is on the order of the distance to the AGN, and these cavities are likely still being or have only recently been inflated by 
the AGN, so the computed rise time is not a reliable estimate of the cavity age. 
s Part of a "split" or "Russian doll" cavity. 

^ Cavity is only marginally detected, tabulated properties may not be reliable. 
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Fig. 1. — Exposure corrected, background subtracted 0.3-2 keV Chandra image of NGC 5813. The image has been smoothed with a 1.5" 
radius Gaussian and point sources have been filled-in by randomly drawing from a Poisson distribution fit to a local background annular 
region. The image shows two pairs of cavities, plus an outer cavity to the northeast, two sharp edges to the northwest and southeast, and 
bright rims around the pair of inner cavities. 
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Fig. 2. — The core from the image in Figure^ unsmoothed. The point-like AGN is visible, offset ^ 0.5 kpc southeast of the axis of 
symmetry of the inner cavities. There is a sharp edge just southeast of the AGN, and the two rings surrounding the inner cavities overlap 
to form an indented structure to the northwest. The black circles define the inner edges of the intermediate cavities. 
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Fig. 3. — Same as Figure [T] but smoothed with a 3" radius Gaussian to show faint features in the outer regions. The image shows a 
weak outer cavity to the southwest, and a faint outer edge-like feature (indicated by the multiple short arrows) that encircles NGC 5813. 
To the northeast the edge structure is less clear, as there is brighter diffuse emission extending across the edge out to larger radii in this 
direction. 
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Fig. 4. — 0.3—2.0 keV Chandra image from Figure [l] (with sources included), divided by a fitted beta model and smoothed with a 6" 
Gaussian to more clearly show surface brightness fluctuations over a wider dynamic range. The cavity regions listed in Table |3] are overlaid 
for clarity. We see the two small inner cavities at ~ 1.5 kpc, the middle cavities (with the northeastern "Russian doll" cavity split into two 
piece) at ~ 8 kpc, and the outermost cavities at ~ 20 kpc. 
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Fig. 5. — Chandra 0.3-2.0 keV image of the core, with 1.36 GHz VLA B array configuration (blue) and 235 MHz GMRT (green) contours 
overlaid. The contours start at 3cr and are spaced by a factor of two. Low frequency radio emission fills the intermediate cavities near the 
edge of the FOV, and emission at both frequencies fills the inner cavities. The central peak of the 1.36 GHz contours is coincident with 
the AGN, while the peak of the 235 MHz contours is northwest of the AGN, roughly coincident with the center point of the elliptical edge 
defined by the inner shock front. The dashed lines indicate the beam sizes for each radio observation. 
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Fig. 6. — Temperature map derived from the ACIS-S3 data, with the cavity regions from Figure |4] overlaid. The color-bar gives the 
temperature in keV. The shocks stand out as temperature peaks coinciding with the surface brightness edges. Cooler gas follows the 
southeastern edges of the intermediate surface brightness cavities indicated in Figure [l] 




Fig. 7. — Pseudo pressure (left) and entropy (right) maps, in arbitrary units. The pressure map was calculated as kTA^^^ and the 
entropy map as kTA~^^^ , where A is the APEC normalization scaled by the area of the extraction region. The pressure jumps at the 10 kpc 
shocks are visible ~ 10 kpc northwest and southeast of the central peak. 




Fig. 8. — A high-spatial resolution temperature map of the core of NGC 5813, with Chandra X-ray 0.3-2.0 keV logarithmic surface 
brightness contours overlaid. The color-bar gives the temperature in keV. The bright rims surrounding the innermost cavities are revealed 
to contain hot gas, probably shock-heated by the central AGN during a recent outburst. 
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Fig. 9. — Pseudo pressure (left) and entropy (right) maps, in arbitrary units, derived as in FigurejT] X-ray surface brightness contours are 
overlaid in black. There is a sharp pressure jump coincident with the temperature jump southeast of the core seen in Figure |8] identifying 
this feature as a shock. Also visible are local pressure minima (and entropy maxima) at the location of the inner cavities, presumably due 
to the additional non-thermal pressure support from the radio emitting plasma in the cavities. 
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Fig. 10. — Azimuthally averaged radial profiles for (top to bottom) temperature, electron density, pressure, and entropy, extracted from 
circular annuli. The top panel shows the projected (black circles) and deprojected (red triangles) temperature profiles. All other panels 
show deprojected values. The vertical dashed lines mark the average radius of the 1.5 kpc (left) and 10 kpc (right) elliptical shock fronts 
indicated in Figure m The shocks are smeared over multiple bins due to their non-spherical structure. The temperature and pressure 
increase at each shock front. The temperature profiles show a central temperature spike, even though the region of the central AGN has 
been excluded from the fits. This is because the profile center is in the region of the hot overlapping rims of the central cavities, shown in 
Figure |2] The dips in pressure and density at ~6 kpc are duo to the X-ray cavities. 
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Fig. 11. — Projected (black circles) and deprojected (red triangles) temperature profiles across the northwestern (left) and southeastern 
(right) shocks. Distance is measured from the center of curvature of the corresponding edge and should not be c ompared with distances in 
the azimuthally averaged profiles. The dashed lines mark the positions of the density jumps calculated in § |5.1[ 
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Fig. 12. — Integrated emission measure profiles across the northwest (left) and southeast (right) edges around the middle cavities, 
indicated in Figure^ Distances are measured from the center of curvature defined by the shock, not from the central AGN. Each profile 
shows a sharp edge at ~13 kpc. The solid lines show the projected best-fitting discontinuous power-law density models, and the vertical 
dashed lines mark the best-fitting break radii corresponding to density discontinuities. 
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Fig. 13. — Integrated emission measure profile across the 1.5 kpc edge just southeast of the core, most clearly visible in Figure[2] which 
is coincident with the edge seen in the core temperature map (Figure [8]|. 
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Fig. 14. — Temperature profile from our fiydrodynamical simulations of a point explosion in an isothermal sphere, chosen at a time when 
the corresponding surface brightness profile best matches the observed profile for the northwestern outer shock. The solid lines show the 
emission-weighted projected temperature profile, with 90% confidence intervals. The points are from results of XSPEC simulations, with the 
model fold ed through the Chandra response and projected. The predicted projected temperature rise is consistent with what is seen in 
Figure [TT] even though the true temperature rise in the model is larger and consistent with the prediction from ID shock jumps conditions 
based on the Mach number. 



28 



RANDALL ET AL. 



10 



-5 




150 200 
R (arcsec) 



Fig. 15. — 0.3-2.0 keV background-subtracted surface brightness profiles in two wedges across the outer edge indicated in Figure [3| to 
the northwest (red circles) and south (blue triangles). The northwestern profile was measured between 12° - 79° and the southern profile 
between 206° - 320° (angles measured from west to north). The northwestern profile shows a jump, followed by a change in slope, at ^ITO" 
(vertical dashed line), roughly the location of the outer edge. The southern profile also shows a change in slope at ^lyO", though with no 
associated jump. 
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Fig. 16. — 0.3-2.0 keV background subtracted surface brightness profiles across the outer cavities indicated in Figure 13] For the 
northeastern cavity (left), the profile was extracted between 112° - 153° (measured from west to north), and between 286° - 310° for the 
southwestern cavity (right). Each profile shows a significant dip at the location of the cavity. For the northeastern cavity, the bright outer 
rim is also visible. 
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Fig. 17. — Left: Ho image of the core of NGC 5813 taken with the SOAR telescope. The Chandra X-ray contours are overlaid in green. 
The H« filaments extend to the southern edges of the intermediate cavities, and are co-spatial with the trail of cool gas seen in the X-ray 
temperature map in Figure |6] In the center, the Ha emission anti-correlates with the inner X-ray cavities. The image shows the cool Ha 
gas being displaced as the central cavities are inflated, and lifted by the intermediate cavities as they rise buoyantly. Right: Close-up of 
the central part of the temperature map shown in Figure |6] with Ha contours overlaid in white. The Ha is co-spatial with the cool X-ray 
filament. 



